Low Temperature Magnetic Properties of the Double Exchange Model 
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We study the ferromagnetic (FM) Kondo lattice model in the strong coupling limit (double exchange 
(DE) model). The DE mechanism proposed by Zener to explain ferromagnetism has unexpected 
properties when there is more than one itinerant electron. We find that, in general, the many-body 
ground state of the DE model is not globally FM ordered (except for special filled-shell cases). Also, 
the low energy excitations of this model are distinct from spin wave excitations in usual Heisenberg 
ferromagnets, which will result in unusual dynamic magnetic properties. 



The double exchange (DE) model JT] ||] has attracted 
much recent attention p]-|l"l| because of its anticipated re- 
lation to the Mn-oxide perovskites La\- x A x MnOz (A — 
Co, Sr, Ba) materials with colossal magnetoresistance 
(CMR) fll2|- jlqj . To explain ferromagnetism in these Mn- 
oxides, Zcncr |lj] introduced a DE mechanism, in which 
local S — 3/2 spins of the three Mn ti g <i-electrons be- 
come ferromagnetically coupled due to coherent hopping 
of e g electrons. 

The DE mechanism can be derived from the follow- 
ing ferromagnetic (FM) Kondo lattice model Hamilto- 
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where a is the Pauli matrix and Si is the local spin S — 
3/2 of three Mn t^g (i-electrons. The operators c ia (cj ) 
annihilate (create) a mobile e g electron with spin a at site 
i. The DE mechanism requires the Hund's rule coupling 
J H > t/S. 

The Hamiltonian (Q) was studied by Anderson and 
Hasegawa H using a single electron on two sites. They 
concluded that, in the low-energy band, the electron spin 
must be parallel to the local spin at every site, leading 
to a renormalization of the electron hopping integral by 
{Sq + \)/ (25* + 1), where Sq is the total spin of the sub- 
system consisting of the two localized spins on sites i and 
j and the electron. Based on this single electron result, 
in most theoretical calculations the original model (Q) is 
mapped onto a non-interacting model of spinless fermions 
with spin disorder scattering . In this non- interacting 
model, the low temperature phase is ferromagnetic (FM), 
due to the gain of electron kinetic energy if the local spins 
are aligned. This DE mechanism is also implemented in 
many-body treatments |i|,||f7|,[l0| by using FM mean-field 
solutions or by expanding around them. However, in 
the original model (||) the ground states are not FM for 
some fillings due to the Pauli principle. This is read- 
ily seen in the half-filling case (see below), in which the 
ground state has antiferromagnetic (AFM) correlations. 



In this Letter, we study the Hamiltonian ([!]) using ex- 
act diagonalization and variational wavefunctions. We 
show that the DE mechanism is a purely single electron 
property, and that the many-body ground states are typ- 
ically not FM JIt]] , except for special fillings in which the 
non-interacting states at the Fermi level are completely 
filled. The ground states in model (0) are generally only 
locally FM without global FM ordering. Also, the low en- 
ergy excitations are different from spin wave excitations 
in Heisenberg FM systems even if the ground state has 
FM ordering at special fillings. 

The Hamiltonian (Q) conserves total momentum K, 
total spin S\ and S l z . For the FM states, S 1 = S l max = 
NS + n/2, where N is the number of sites, and n is the 
number of electrons. It can be easily seen that the wave- 
functions of low energy FM eigenstates are simply the 
Slater determinants of spinless fermions \k±, k,2, ■ ■ ■ , fen)® 
\FM), with energy J2i(~JhS + eg ), where is the 
non-interacting electron dispersion. \FM) denotes the 
FM states for local spins, with the electronic spins be- 
ing fixed by local spins. We will denote the lowest FM 
eigenenergy by Eq. For our numerical calculations, we set 
t = 1, Jh = 20. In the manganese perovskites S = 3/2, 
however we will take S = 1/2 for most of our numerical 
calculations; some calculations are repeated for smaller 
systems at S = 3/2, and we obtain qualitatively the same 
results. 

For one electron (n — 1), it can be shown analyti- 
cally that in the ground state 5* = 5'^ laa , for Jh > 0, 
which is consistent with the DE mechanism. At half- 
filling n = N, there is an induced AFM coupling ~ t 2 / Jh 
between spins due to the Pauli principle: the electrons 
cannot hop in the FM state; they can only hop if the 
neighboring sites are empty or have anti-parallel spin, 
inducing an AFM coupling. Thus the ground state at 
half-filling is AFM jL8|. The question now is, what hap- 
pens for general fillings. We can take some guidance from 
a heuristic non-interacting electron picture: the electrons 
have to be parallel to the local spin at every site, since 
JnS/t S> 1. Due to the renormalization of the elec- 
tron hopping integral t — > t{S>Q + h)/(2S + 1), mobile 
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electrons have the largest bandwidth in the FM states; 
electronic spins in these states are totally polarized, and 
then fill states up to the Fermi level. Non-FM states of 
the local spin will reduce the bandwidth; but the elec- 
tronic spins are not totally polarized, and so can have a 
lower kinetic energy. Thus it is plausible that at some 
fillings the ground state can be non-FM. It is straightfor- 
ward to see that non-ferromagnetic ground states exist 
in the limit where the local spins are classical (S — > oo). 
Consider a ID ring where the local spins uniformly cir- 
cle the north pole at an angle 9. The magnitude of the 
effective hopping t e ff is reduced as 9 increases, but an 
additional magnetic flux (Berry phase) is simultaneously 
added. It is easy to verify that for large Jh the change 
in energy AE = —cit 2 sin 2 (9)/ Jh, where c\ is a positive 
constant for an even number of electrons. For finite Jh, 
the Berry phase reduces the total energy by more than 
t e f f increases it, and the spins prefer to lie on the equator 
& 

That this mechanism is also relevant for small S can be 
shown for two electrons in any dimension, by construct- 
ing a wavefunction that is not FM but lower in energy 
than Eq. We denote the zone center fco = with (lowest 
Bloch) energy eo, and the first excited states with energy 
ei by Qi, with i = 1, • • • , M (M > 2). Then we can write 



a wavefunction with S = SI 



1 as follows: 
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where |0) is the electronic vacuum state. The variational 
energy E v decreases with increasing Jh ■ As Jh /t — * oo, 
x -> 1, and E v = E - (ei - e )(M - 1)/(2SN + M + 
1). If the number of sites N > 2, the energy difference 
Eq — E v is finite as Jh — > oo. Our exact diagonalization 
calculation shows that the ground state is S l — in ID 
and 2D for two electrons (n = 2). Consequently, the 
trial wavefunctions can only tell us that there exist lower 
energy non-FM states. These trial wavefunctions have 
finite probability of two electrons being in |fco = 0) by 
using the degeneracy of the non-interacting states at the 
Fermi level. 

In ID, we can extend the results for n = 2 to general 
fillings: The ground states for even numbers of electrons 
are non-FM for any n — 2m (m integer) and N. This 
can be proved by constructing the following wavefunction 
with S** = S* - 1 Poll: 



l*> = y er lqal {S7 + xct c n ) \vac) - — * - 
I / £_j v i n ipi i 2NS + 2m 



with | vac) = c[ 1 _ m)9oT ---cJ mftT |0)|FAf) and q = 
±2ir/N. As for n = 2, the variational energy E v de- 
creases with increasing Jh- It can be shown that, for 



Jh /t — > oo, x — ► 1, and E v < Eq- So there is at least one 
eigenstate with S l = S^^ — l, K = (m—l)q with energy 
less than E . We believe that the difference between even 
and odd number of electrons arises as follows: In the FM 
states, the electrons fill the non-interacting states up to 
the Fermi level; for an odd number of electrons, the ±fc/ 
states at the Fermi level are all filled; for an even num- 
ber of electrons, only one of the ±fc/ states is filled. This 
kind of degeneracy is related to the possibility of having 
non-FM ground states. The wavefunctions (|2|) and (|^) 
have a lower energy than Eq due to this degeneracy. In 
ID systems, we calculated the ground state and some low 
energy excited states using Lanczos methods for sizes up 
to N = 16, n = 2; and N = 10, n = 6. We find that the 
ground state is always S l = and K = it for even num- 
bers of electrons, and FM and K = for odd numbers of 
electrons. 
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FIG. 2. Site indices for the 2D 4 x 4 lattice. 

The trial wavefunction in ID cannot be generalized to 
higher dimensions except for n — 2. However, we believe 
that the origin of the non-FM ground state should be the 
same for any dimension. Recall that our results in both 
the n = N and the n = 2 cases are independent of dimen- 
sionality. We speculate that the ground state of model 
(|l|) is non-FM in higher dimensions whenever the single- 
particle states at the Fermi level of the FM ground state 
are not completely filled. This has been confirmed in our 
finite size exact diagonalization calculations. The filling 
of the non-interacting single electron states plays an im- 
portant role here because the wavefunction of the FM 
ground state is just a Slater determinant wavefunction of 
the spin polarized non-interacting electron system. From 
the single electron dispersion of the non-interacting sys- 
tems (Fig. (|l|)), we see that the filled-shell cases are 1, 7 
for the \/8 x V8 lattice and 1, 5, 9 for the VW x VlO 
lattice. For these fillings, we find that the ground states 
are FM. The non-fillcd-shcll ground states have 
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S = (2), 1/2 (3), (4), 1/2 (5), (6), (8), for N=8, 
5 = (2), 5/2 (3), 4 (4), 7 (6), for N=10 and 
S = 0(2) for N=16. 

(The number in parenthesis is the number of electrons). 
In all these cases the ground states are non-FM. We also 
calculated the 5 = 3/2 case up to N = 8 and n = 4. The 
results are similar to S = 1/2. One different feature in 
a 2D (e.g. y/lO x vlO) lattice is that the ground states 
are not necessarily at the smallest 5* for non-filled shell 
fillings. But for all the cases we have checked, there is 
at least one state in each 5* < which has a lower 

energy than E . 

What is the magnetic structure of the non-FM ground 
state? At n = TV, the ground state is AFM as described 
above. At n = 2, the ground state also has S 1 ' = 0, but 
the spins have local FM ordering, which is evident in the 
correlation functions (see Table | for a 4 x 4 2D lattice) . 
For other fillings with non-FM ground states S* <C ^maxi 
the magnetic structures are most likely to be locally but 
not globally FM ordered (see Fig.(||)). For ID systems, 
the local spin structures are spiral with a pitch equal to 
the length of the chain. This can be seen in the correla- 
tion function ((SixS2)-(SiXSi+i)), which is only weakly 
dependent on i and has magnitude ~ (0.25 sin(27r/A r )) 2 . 
Preliminary numerical results suggest that the 2D non- 
FM ground states are also noncolinear for non-filled-shell 
cases. The details of the spin structures for ID and higher 
D are not the major concern here and will be discussed 
elsewhere. What we can conclude from these correlation 
functions (see Fig.(^) and Table |) is that the spin-up 
electrons are repelled from spin-down electrons and the 
local spins with greatest separation are anti-parallel: thus 
there is local (i.e. FM domains), but not global FM or- 
dering. 



wave" (SW) excitations. However, there are many other 
branches of AS* = 1 excitations in model ([!]). For the 
quasi SW states, the dispersion is not a cosine, as it 
would be in a near-neighbor Heisenberg ferromagnet (see 
Fig.(||)). In Fig.(^) we show the dispersion of SW states 
of a ID system with N = 16, n = 1,3,5. From Fig.(g), 
we can see that the effective spin-spin coupling J softens 
at short wavelength, and the softening decreases with in- 
creasing electron density. This softening at short wave- 
length (large wavenumber) is due to a spin polaron effect: 
the electron density is not uniform and is smallest at the 
site with a flipped-spin. For long wavelength SW excita- 
tions, the electron density is nearly constant, so the SW 
dispersion is quadratic. As K increases, the electrons in- 
creasingly avoid the flipped spin, so the effective spin-spin 
coupling J becomes smaller. One can show using a Bcthc 
ansatz that the n = 1 excitations shown in Fig.(|J) cannot 
have an energy greater than 2£[1 — cos(7r/iV)]. This is \ of 
the lowest electronic excitation energy, which promotes 
a free electron from K — to K — 2-k/N . 
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FIG. 4. Dispersion of SW states at TV = 16, n = 1,3,5. 
We have rescaled the n = 1, 3 curves. The solid line is the 
cosine curve fitted by the lowest two points. Parameters: 
S = 1/2, t = 1. J H = 20. 
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FIG. 3. Correlation functions (nln-"-^, (TiinJ + i), 
(Si ■ S,-+i) for the ID N = 10, n = 4, S = 1/2 system. 
The curves are rescaled to fit on the same graph, so only the 
relative magnitudes are representative. 

The excitation spectrum of model ([!]) is also un- 
usual. For the filled-shell cases, there are quasi "spin 



FIG. 5. Dispersion of low energy states at N = 8, n = 2. 
Symbols: filled O: S = 0; •: S = 1; O: S = 2; O S = 3; □: 
S = 4; A: S = 5. Parameters: S = 1/2, i = 1, J H = 20. 

For the non-filled-shell fillings, the low energy excita- 
tion spectrum is complicated (see Fig.(||)). As discussed 
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above, the ID spin structure in non-fhled-shell cases has 
spiral ordering. This might suggest that there will also 
exist quasi SW low energy excitations. However, from 
Fig.([|), we can see that it is difficult to interpret the 
excitations as spin waves. 

We have shown that model (|l|) exhibits non-FM ground 
states on finite lattices for non-fillcd-shcll cases in both 
ID and 2D. (This has also been confirmed in 3D for up to 
4 electrons on a 3 x 3 x 3 lattice.) There are many small 
S f states with lower energy than the FM ground state 
energy E (see Fig. (El)). It is believed that there is weak 
AFM super-exchange coupling between local spins in the 
CMR materials 0Jq,^ll . The low temperature phase will 
further deviate from FM if this is included. Experiments 
on CMR materials show that the low temperature phase 
is FM at dopings ~ 0.2 — 0.4 with a quadratic long wave- 
length spin wave dispersion at low temperatures |p2| , p3f . 
One possibility is that the model Hamiltonian (JjJ) ioes 
not describe the physics in the CMR materials: a more 
complete model should include, e.g. two interacting or- 
bitals and electron-phonon (e-p) couplings. Another pos- 
sibility is that the locally FM ordered phase found here is 
experimentally indistinguishable from global FM order- 
ing. Because of the inhomogeneity of electronic density 
due to disorder or e-p coupling, the filled-shell phases co- 
exist with non- filled-shell phases. Thus there are quasi 
SW excitations accessible to neutron experiments. This 
is of course speculative, but some of the unusual exper- 
imental magnetic properties [^3] in CMR materials 
might be explained along this line. Notice that the SW 
excitations of model (Q) in filled-shell cases is only one of 
the many branches of AS — 1 excitations; other branches 
can have considerable higher energies. Thus T c in model 
(|l|) may be higher than estimated using the spin stiff- 
ness from SW excitations alone, i.e. the spin coupling is 
energy dependent. 

In a homogeneous system, the local FM spiral ordering 
discussed above might be experimentally indistinguish- 
able from the global FM ordering. However, if a double 
exchange system is inhomogeneous and can be separated 
into small mesoscopic systems due to Anderson localiza- 
tion, the spiral related states will have distinct properties 
from the FM states. To verify that non-FM states occur 
not only for periodic boundary conditions, we have cal- 
culated the ground states for various small clusters using 
open boundary conditions. We find that the ground state 
can still be non-FM in D > 1 finite clusters (2x3,3x4, 
etc) with open boundary conditions. 

In this work, we showed that the DE mechanism is 
a single particle mechanism which does not give FM in 
many-electron systems. In the DE model (Q), there are 
many small S* states with lower energy than the FM 
ground state energy Eq, except for filled-shell cases. In 
ID, for an even number of electrons, the ground state 
has S l = and momentum K = tt. In higher D, the 
ground state is FM only if the non-interacting states at 



the Fermi level are completely filled. Otherwise there are 
states in each S" < S^, aa . subspace, with energy < E a . 
The spin correlations show that the non-FM low energy 
states have local FM-domain-like structures (spiral order- 
ing in ID): two groups of electrons (spin-up, spin-down) 
tend to avoid each other by forming two ferromagnetic 
domains. We suggest that this low temperature phase 
might be relevant to the CMR materials, which have been 
shown to possess unusual magnetic properties in recent 



experiments 23 - 25 1 . 
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TABLE I. The ground state static correlations functions 
for 2D N = 16 (4 x 4), n = 2, S = 1/2. {Si ■ Sj) is the 
correlation function of the total spin (local + itinerant). Site 
indices are shown in Fig. (2). 
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